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ABSTRACT 

Many important techniques for investigating the properties of extragalactic 
radio sources, such as spectral-index and rotation-measure mapping, involve 
the comparison of images at two or more frequencies. In the case of radio inter- 
ferometric data, this can be done by comparing the CLEAN maps obtained at 
the different frequencies. However, intrinsic differences in images due to the 
frequency dependence of the radio emission can be distorted by additional 
differences that arise due to source variability (if the data to be compared is 
obtained at different times), image misalignment, and the frequency depen- 
dence of the sensitivity to weak emission and the angular resolution provided 
by the observations (the resolution of an interferometer depends on the lengths 
of its baselines in units of the observing wavelength). These effects must be 
corrected for as best as possible before multi-frequency data comparison tech- 
niques can be applied. We consider the origins for the afore-mentioned factors, 
outline the standard techniques used to overcome these difficulties, and de- 
scribe in detail a technique developed by us, based on the cross-correlation 
technique widely used in other fields, to correct for misalignments between 
maps at different frequencies. 



1 INTRODUCTION 

The radio emission of core-dominated, radio-loud Active Galactic Nuclei (AGN) is syn- 
chrotron radiation generated in the relativistic jets that emerge from the nucleus of the 
galaxy, presumably along the rotational axis of a central supermassive black hole. One im- 
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portant source of information about the physical conditions in the radio-emitting regions is 
the distribution of the spectral index a over the source {Si, cx i^", where Si, is the source 
flux at frequency u). The core region is typically observed to be at least partially optically 
thick, with a nearly flat or inverted spectrum, while the jets are optically thin, with nega- 
tive spectral indices. The spectrum may also flatten in regions of the jet in which there is 
re-acceleration of electrons or low-frequency absorption (e.g. Gabuzda, Pushkarev & Gar- 
nich 2000; Gabuzda, Gomez & Agudo 2001). Synchrotron radiation can be highly linearly 
polarised, to ~ 75% in the case of a uniform magnetic field (Pacholczyk 1970), and linear 
polarisation observations can yield unique information about the orientation and degree of 
order of the magnetic field in the synchrotron source, as well as the distribution of thermal 
electrons and the magnetic- field geometry in the immediate vicinity of the AGN (e.g., via 
Faraday rotation of the plane of polarisation) . 

The compact radio emission of such AGN can be probed with high resolution using Very 
Long Baseline Interferometry ( VLBI) . The radio telescopes in the interferometric array can 
be separated by hundreds or thousands of kilometres, making it infeasible to physically link 
(synchronise) them electronically, and high-accuracy timing signals must be recorded to- 
gether with the data, so that the signals obtained at different antennas can be accurately 
synchronised during correlation. In practice, the amplitudes and, especially, phases of the 
measured complex visibility data unavoidably contain unknown errors, which can conve- 
niently be expressed via antenna-based complex gain factors: 

where V^^^ and V^J'^'^ are the observed and true visibility functions on the baseline between 
antennas i and j, and Gj, Gj are the complex gain factors for antennas i and j. The complex 
gains G must be determined and removed from the data in order to achieve the highest- 
quality images possible for the radio-telescope array used. This is normally done iteratively, 

via alternate application of self-calibration (Fort and Yee 1976, Cotton 1979, Readhead & 
Wilkinson 1978, Readhead et al. 1980, Cornwell & Wilkinson 1981) and a deconvolution 
method, such as CLEAN (Hogbom 1974). 
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2 MATCHING AND ALIGNING IMAGES AT DIFFERENT 
FREQUENCIES 

Since AGN are variable, it goes without saying that multi-frequency data to be compared 
must correspond to epochs separated by time intervals appreciably less than the timescales 
for variability of the source. When preparing data for techniques involving comparison of 
multi-frequency data, various instrumental differences between the datasets must also be 
taken into account. VLBI datasets at different frequencies will have different angular res- 
olutions and sensitivities to structures on various scales due to the different baseline cov- 
erages of the observations. One approach to reducing these differences when comparing 
multi-frequency images is to match the baseline coverages at different frequencies by giving 
relatively low weights to the longest baselines at higher frequencies and to the shortest base- 
lines at lower frequencies, e.g., via tapering of the visibility data. Alternatively, the images 
can be obtained without such weighting, but then all be convolved with the same CLEAN 
beam before comparison. The size of the CLEAN beam to be used in this case is ordinarily 
roughly equal to the central lobe of the dirty beam for the lowest frequency. The accuracy 
to which the positions of the CLEAN components arc known is limited by the resolution of 
the observing system, and convolving with a beam that was much smaller than the central 
lobe of the dirty beam would result in a "superresolved" image that may not be reliable. 

Iterative imaging via self-calibration and a deconvolution algorithm such as CLEAN is 
generally quite effective, but the absolute information about the coordinates of the source 
on the sky is lost during phase self-calibration, which essentially places the centre of gravity 
of the radio brightness distribution at the phase centre, which has coordinates (0,0). Be- 
cause most radio-loud AGN are highly core-dominated, directly comparing multi-frequency 
self-cahbrated VLBI images of AGN essentially amounts to aligning these images on the ob- 
served VLBI core, which usually coincides very closely with the peak of the radio brightness 
distribution. 

However, the standard theory of extragalactic radio sources (e.g. Blandford & Konigl 
1979) predicts a frequency-dependent shift in the location of the VLBI core due to opacity 
effects in the core region. Reabsorption of synchrotron radiation takes place in the ultra- 
compact region near the central engine of an AGN, a mechanism which is more efficient 
at low frequencies. Consequently, the peak brightness appears further along the jet axis 
in lower frequency observations. Thus the alignment of multi-frequency images on their 
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VLBI core results in a misalignment between images of different observing frequencies. This 
prediction is supported by observation: the frequency-dependent shift in core position has 
been measured for several quasars and micro-quasars (see Lobanov 1998 and references 
therein), and the effect is discussed in terms of its impact on high precision astrometry by 
Fey (2000), Chariot (2002), Ros (2005) and Boboltz (2006). Lobanov (1998) explains in 
detail how the frequency-dependence of the shift depends on physical conditions near the 
central engine. 

It is thus necessary to correctly align images prior to applying multi-frequency data- 
analysis techniques. This can be achieved in one of two ways. The first is by phase-referenced 
observations, first employed by Marcaide & Shapiro (1984), in which a nearby source (or 
sources) is observed along with the target source. The reference source would ideally be 
a point source, to eliminate structure effects including those discussed above, but this is 
rarely possible since most sources show extended structure on the milliarcsecond scales 
available with VLBI. The position of the target source relative to the reference source can 
then be determined. The second method involves aligning images according to the positions 
of optically thin jet components (i.e. components optically thin to synchrotron radiation, 
so that their positions are not affected by absorption effects such as those occuring in the 
core) that are present in both images to be compared. This can be non-trivial, particularly 
if the source has a complicated structure, but has been employed effectively by several 
authors, such as Paragi et al. (2000), who used this method to determine the radio core 
shift in 1823-1-568. This difficulty in aligning complex images without distinct optically thin 
components detected at all frequencies to be compared was the main motivation for us to 
consider alternative methods of image alignment. 



3 IMAGE ALIGNMENT VIA CROSS-CORRELATION 

The method we have developed to align multi-frequency images is based on the cross- 
correlation technique widely used in many fields, including biomedical signal processing and 
imaging (Panescu 1993, Prank & McEwan 1992) and remote sensing (Hartl 1976). Cross- 
correlation provides a measure of how closely correlated (i.e. how alike) two functions are. 

The use of this measure in image alignment gives an objective, quantitative assessment of 
how well two images are aligned, and does not depend on the presence of very compact 
features. By applying different shifts between images and calculating the cross-correlation 
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coefficient for each shift, it is possible to determine which shift results in the best alignment. 
The normalised cross-correlation coefficient we used (see, e.g., Dunn & Clark 1974) is defined 
in two-dimensions as: 



where n is the number of pixels in each direction in the two-dimensional images to be 
compared, Iui,ij and Iu2,ij ai'c the intensities for the maps at frequencies ul and z/2 at pixel 
{RAi, Decj), and I^i and 1,^2 are the mean values of these two intensities over the region 
analyzed. Although the source emission at a given location varies with frequency, all the 
radiation observed at radio frequencies is believed to be synchrotron radiation from the same 
population of relativistic electrons. Therefore the optically thin total-intensity (/) structures 
should be well-correlated despite local changes in the spectrum of the synchrotron emission 
due to variations in the local magnetic field, interaction with the surrounding medium and 
other effects. Using this (reasonable) assumption, it is possible to determine the shift to be 
applied between maps by comparing the structures of optically thin regions of the source at 
different frequencies. The highest correlation between dual-frequency images should therefore 
be obtained when the areas being compared correspond to the same physical region of the 
sky. This method has the advantage that it takes account of the optically thin emission from 
the entire source, not just isolated compact components. 

4 IMPLEMENTATION OF THE CROSS-CORRELATION TECHNIQUE 

The most widely used software for the calibration, imaging and analysis of radio interfero- 
metric data is the National Radio Astronomy Observatory AIPS (Astronomical Image Pro- 
cessing System) package. We have written a C program to implement the cross-correlation 
technique, which is external to but compatible with the NRAO AIPS package. The input 
to the program are two images in the format produced by the AIPS task IMTXT, and 
the program outputs files which can be imported back into AIPS using the task FETCH. 
IMTXT allows the user to export an AIPS image as a text file containing an array of floats 
representing the map values at each pixel location. Conversely, FETCH imports a text file 
in the format exported by IMTXT as a map ffie which can be displayed by any AIPS tasks 
that work with images (e.g. KNTR, TVLOD). 

The input images to the program must have the same pixel size and numbers of pixels 
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Figure 1. Implementation of image alignment by cross-correlation. A sub-area (shaded region) of the first map (outlined in 
bold) is compared with the overlying region of the second map. This sub-area remains constant, while the corresponding region 
of the second map changes as it is shifted relative to the first map (some possible shifted positions are outlined with dashed 
lines). The cross-correlation coefficient between the two regions, given by equation (f), is computed each time to provide a 
measure of how well the areas are correlated. The method assumes that the highest correlation is achieved when the areas being 
compared refer to the same physical region of the sky. The circle labeled "core" represents a region with a specified elliptical 
or circular shape (corresponding to the beam shape) coincident with the position of the optically thick core, that is omitted 
from Map 1 during the cross-correlation calculations. 

(i.e. so that both the size of a pixel and the overall size of the images correspond to the 
same area on the sky), convolved with the same beam, and exported from AIPS using the 
task IMTXT. The user specifies the maximum trial shift to be applied between the images. 
The procedure used to calculate the shift that best aligns the input images is as follows. 

(i) Strips of width A, where A is the maximum shift to be applied to the maps, are 
subtracted from each edge of the first map (Map 1). An area of the same shape as the 
restoring beam, but whose size is specified by the user, is removed from the area in Map 1 
surrounding the (at least partially) optically thick core, whose position will be frequency- 
dependent and which usually corresponds to the peak of the map. In comparatively unusual 
cases, the brightest feature may not correspond to the core; to enable the program to produce 
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reliable results for this situation, we have included the option of the user specifying the 
position of the core. 

(ii) The remaining area is used in the comparison (see Fig. 1). The second map (Map 2) 
is shifted so that different regions overlay this area each time. Since the area of Map 1 is not 
changed, the normalised cross-correlation coefficients can be compared directly to determine 
which part of Map 2 is best correlated with the selected region. 

(iii) Map 2 is first shifted so that its bottom left hand corner corresponds to the bottom 
left hand corner of the selected sub- area of Map 1. This corresponds to the maximum 
negative shift (—A, — A) applied. Map 2 is then shifted in right ascension, one pixel at a 
time, and the cross-correlation coefficient computed each time, until the maximum positive 
shift A is reached. The image is then shifted by one pixel in declination, and the correlation 
coefficients for the next row computed in the same way. This is repeated until the maximum 
shift in both directions (A, A) is reached, which occurs when the top right hand corner of 
Map 2 corresponds to the top right hand corner of the selected subarea of Map 1. 

(iv) The program outputs to the screen the maximum correlation Vmax and the shift 
{dRA,dDec) at which it occurs. 

(v) Having found the best initial alignment between the two images, the program now 
applies the corresponding shift to Map 2, constructs a spectral- index map and blanks any 
optically thick points in this map, taken to be points with spectral indices a > 0. The 
theoretical limiting optically thick spectral index o; = | is rarely reached, but it is generally 
accepted that a positive spectral index implies some optically thick contribution in that 
region. The optically thick regions are blanked after first calculating the shift with only the 
core region blanked, since the shift between the maps can result in a significant change in the 
spectral-index distribution, and thus in those regions which have a > 0. The cross-correlation 
procedure is then repeated, taking into account only the optically thin (i.e. unblanked) 
regions of the source, since the position of these regions is not affected by absorption effects. 

(vi) The program again outputs to the screen the maximum correlation r^ax and the 
shift {dRA,dDec) at which it occurs. Positive shifts correspond to moving the second image 
downward (to the South) and to the left (East) relative to the first image. In other words, a 
feature whose pixel location is {RA,Dec) in the first map is located at {RA+dRA,Dec+dDec) 
in the second map. 

(vii) Three files are output. One is simply a text file containing the array of cross- 
correlation values computed, while the other two are text files of the format recognised 
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by the AIPS task FETCH. One of these hkewise contains the cross-correlation coefficients; 
importing this image into AIPS and plotting it (e.g with KNTR) shows the shape of the 
cross-correlation function. The other image text file displays the subarea of Map 1 that was 
compared; displaying this image can be useful in verifying the area to be blanked around 
the optically thick core. Ideally this area should be large enough to cover virtually all of the 
core region, but not features in the inner jet; it should cover at least one beam size, and 
usually more, since the shape of the beam and the very high flux emanating from the core 
will dominate the VLBI / structure here. 

In all cases we have tested, we found the cross-correlation function to fall off mono- 
tonically from its peak. This lends credence to the hypothesis that the optically thin jet 
structures in the maps should be well correlated when aligned properly (if small-scale vari- 
ations were important the function could show secondary peaks where individual features 
align well rather than falling off uniformly, but we have never found this to be the case). 

Obviously, the accuracy of the estimated shift between a pair of images will depend 
to some extent on the pixel size in the images being compared. In practice, it may be 
expedient, within reason, to use images with slightly smaller pixel sizes in order to derive 
more refined shift estimates. The results should not depend on whether the higher-frequency 
or lower-frequency images is used as Map 1: in either case, it is the relative shift that is being 
determined, and the shifts obtained with the higher-frequency map as Map 1 will be the 
negative of those obtained with the lower- frequency map as Map 1. In practice, it may be 
expedient to compare results obtained blanking out core areas of several different sizes (e.g. 
2.0 times the beam area, 2.5 times the beam area etc.), to ensure that the optically thick 
core region is fully blanked, while as much of the optically thin inner jet is included in the 
correlation analysis as possible. 

5 TESTING/ILLUSTRATION OF THE CROSS-CORRELATION METHOD 

We present two examples of applying the program in practice to align total-intensity images 
of AGN obtained with the NRAO Very Long Baseline Array at 7.9 GHz and 5.1 GHz: (i) 
for 2007-1-777, which displays a fairly distinct optically thin jet feature that could be used to 
align the two maps in the "traditional" way, and (ii) for 2200+420, which displays only fairly 
amorphous optically thin jet emission. The maps compared have matched cell sizes, map 
sizes and beam parameters. The VLBA observations were made over 24 hours in July 2006 
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Superposed 5,9cm & 3.8cm I mops for 2007 + 777: Before olignmenf Superposed 5.9cm & 3.8cm I mops for 2007+777: After alignment 
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Figure 2. Maps of 2007+777 before (left) and after (right) alignment, showing the I contours at 5.1 GHz in green, with / 
contours at 7.9 GHz superposed in black. The convolving beam was 1.80 X 1.55 mas in position angle 87.0°, and was the 
same in all cases. The bottom contour levels are ±0.20, and the contours increase in steps of 1.98. The peak brightnesses 
are 657 mjy/beam (7.9 GHz) and 547 mjy/beam (5.1 GHz). The algorithm has clearly led to a very good alignment for the 
distinct, optically thin feature ~ 7 mas from the core. 

in a snapshot mode, so that the basehne coverage obtained was spread out over the time 
the sources were visible with all or nearly all of the VLBA antennas. The data was obtained 
simultaneously at the different frequencies. The preliminary calibration and imaging of these 
data were carried out in AIPS using standard techniques; some initial results are presented 
by O' Sullivan & Gabuzda (2008), and a more complete analysis is in preparation. 

For each of these two sources, we show a superposition of the 7.9 GHz and 5.1 GHz / 
maps before and after applying the derived alignment shift (Figs. 2 and 4), and the spectral- 
index maps obtained before and after applying the derived alignment shift (Figs. 3 and 5). 
In both cases, we used a cell size of 0.1 mas during the initial total- intensity mapping, but 
used final maps made with a cell size of 0.05 mas as input to the program, to improve slightly 
the accuracy of the relative shifts obtained. 

In the case of 2007+777, the derived shift of the 5.1 GHz relative to the 7.9 GHz image 
was 4 pixels (0.20 mas) to the West and pixels to the South, in the direction of the VLBI 
jet, as expected. The correctness of this shift is immediately obvious via a visual inspection 
of the superposed / images (Fig. 2) and the spectral-index maps before and after applying 
this shift (Fig. 3); the algorithm has obviously aligned a distinct optically thin jet component 
located ^ 7 mas to the West of the map peak. In the case of 2200+420, the calculated shift of 
the 5.1 GHz relative to the 7.9 GHz image was 4 pixels (0.20 mas) to the South and pixels 
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Figure 3. Colour spectral-index maps of 2007+777 before (top) and after (bottom) alignment, with the I contours at 7.9 GHz 
from Fig. 2 superposed. The convolving beam is shown in an upper corner of each image. The shown ranges of spectral indices 
are from —1.89 to +1.35 (top) and from —1.38 to +1.96 (bottom). Spectral artefacts due to misalignment are clearly visible 
in the top spectral-index map (false optically thin emission to the East of the / peak and false optically thick emission at the 
Western end of the jet), which are absent from the bottom spectral-index map. 



to the East, again in the direction of the VLBI jet, as expected. It is less straightforward to 
estimate the correctness of this shift directly from the superposed / maps in Fig. 4, but the 
spectral-index map after applying this shift shows appreciably more regular behaviour, with 
a smooth gradient in the spectral index from North of the core region to the jet extending 
nearly directly to the South; in particular, a spurious region of optically thin emission to 
the North of the peak has disappeared. 

We also show the results of applying the cross-correlation technique to 4.8 and 1.6- 
GHz images of the AGN 1803+784 (Gabuzda & Chernetskii 2003), as an example of the 
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operation of the program when apphed to two images at more widely separated frequencies. 
A superposition of the 4.8 GHz and 1.6 GHz / maps before and after applying the derived 
alignment shift is shown in Fig. 6, and the spectral-index maps obtained before and after 
applying the derived alignment shift in Fig. 7. A cell size of 0.50 mas was used during the 
initial total-intensity mapping, but final maps with a cell size of 0.25 mas were used as input 
to the shift program. The derived shift of the 1.6 GHz relative to the 4.8 GHz image was 5 
pixels (1.25 mas!) to the West and 1 pixel to the North, in the direction of the VLBI jet. 
Although examination of the superposed / images (Fig. 6) does not enable an unambiguous 
estimate of the needed shift "by eye," the correctness of the shift derived by our program 
is immediately obvious via a visual inspection of the spectral-index maps before and after 
applying this shift (Fig. 7). At first glance, the "slanting" boundary between the regions of 
optically thick and thin emission near the core seems suspicious, but in fact, the 4.8-GHz 
VSOP space- VLBI image for this epoch shows that the VLBI jet initially emerges to the 
Northwest (Gabuzda 1999), and this feature likely real (the implied spectral-index gradient 
is roughly perpendicular to the direction of the small-scale jet). 

Finally, we show the two-dimensional plots of the cross-correlation functions output by 
the program for 2007+777, 2200+420 and 1803+784 (Fig. 8). The cross-correlation plots 
contain a single peak, with a monotonic fall-off in the correlation coefficient with distance 
from the derived optimal alignment shift. 

6 CONCLUSION 

We have developed a C program to determine the shift between two VLBI images based on a 
cross-correlation analysis of the images. In the past, it has been necessary to determine such 
shifts for images of AGN by aligning compact optically thin jet components, which can be 
a somewhat subjective and not entirely unambiguous procedure. In addition, this method 
is difficult to apply to complex and/or extended AGN jets without compact optically thin 
features suitable for such an analysis. The great advantage of our new approach is that it 
provides a straightforward, objective means to determine the shift between two images that 
makes use of all optically thin regions in the source structure, not just individual chosen 
features. Our tests have shown that the program produces reliable shifts for images both 
with and without distinct optically thin features. The code can be compiled using a standard 
C compiler, and has been designed to take input files written by the AIPS task IMTXT, 
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Figure 4. Maps of 2200+420 before (left) and after (right) alignment, showing / contours at 5.1 GHz in green, with / contours 
at 7.9 GHz superposed in black. The convolving beam was 2.25 X 1.59 mas in position angle —13.8° , and was the same in all cases. 
The bottom contour levels are ±0.20, and the contours increase in steps of 1.98. The peak brightnesses are 1940 mjy/beam 
(7.9 GHz) and 1520 mJy/beam (5.1 GHz). In this case, there is no compact, distinct optically thin feature on which to base 
the derived shift between the two images, but the algoritm has nevertheless properly aligned the images, as is clear from a 
comparison of the corresponding spectral-index maps in Fig. 5. 

and to output files that can be read by the AIPS task FETCH, making it straightforward 
for radio astronomers familiar with AIPS to implement the code. The code can be obtained 
by contacting D. C. Gabuzda (gabuzda@phys.ucc.ie). 
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Figure 6. Maps of 1803+784 before (left) and after (right) alignment, showing / contours at 1.6 GHz in green, with / contours 
at 4.8 GHz superposed in black as an example of the operation of the cross-correlation technique applied to images separated 
by a larger frequency difference. The convolving beam for both images was 5.49 X 4.49 mas in position angle 51.2°. The bottom 
contour levels are ±0.15 (4.8 GHz) and ±0.28 mjy (1.6 GHz), with the contours increasing in steps of 2.0 in both cases. 
The peak brightnesses are 2270 mjy/beam (4.8 GHz) and 1640 mjy/beam (1.6 GHz). The 4.8 and 1.6 GHz "peaks" of the 
somewhat diffuse jet feature roughly 25 mas from the core are at appreciably different positions in the aligned maps, but an 
inspection of the spectral-index maps in Fig. 7 demonstrates that the algoritmm has properly aligned the overall optically thin 
jet structure. (Maps adapted from Gabuzda & Chernetskii 2003.) 



Gabuzda D. C, Pushkarev A. B. k Garnich N. N., 2001, MNRAS, 327, 1 

Hartl P. 1976, in Remote Sensing for Environment Sciences, Ed. E. Schanda (Springer- 

Verlag), 304 
Hogbom J. A. 1974, A&AS, 15, 417 
Lobanov A. P. 1998, A&A, 330, 79 
Marcaide J. M. & Shapiro I. I. 1984, ApJ, 276, 56 

O'SuUivan S. P. & Gabuzda D. C. 2008, in Extragalactic Jets: Theory and Observation 
from Radio to Gamma Ray, Ed. T. A. Rector & D. S. De Young (Astronomy Society of 
the Pacific), in press. 

Pacholczyk A. G. 1970, Radio Astrophysics (San Francisco: Freeman and Co.) 

Panescu D. 1993, in Biomedical Digital Signal Processing, Ed. W. Tompkins (Prentice Hall), 
216 

Paragi Z., Fejes I. & Frey S. 2000, in International VLBI Service for Geodesy and Astrome- 
try 2000, General Meeting Proceedings, Ed. N. R. Vandenburg & K. D. Baver, NASA/CP- 
2000-209893, p. 342 



Aligning VLBI Images of AGN 15 




10 -10 -20 -30 

MilliARC SEC 

Grey scale flux ranges -1.500 0.500 SP INDEX 
Cont peak flux = 1.6403E4-00 JY/BEAM 
Levs = 1.640E-02 * (-0.280, 0.280, 0.550, 1.120, 
2.240, 4.480, 8.960, 17.92, 35.84, 71.68) 




10 -10 -20 -30 

MilliARC SEC 

Grey scale flux range= -1.500 0.500 SP INDEX 
Cont peak flux = 1.6403E+00 JY/BEAIVI 
Levs = 1.640E-02 * (-0.280, 0.280, 0.560, 1.120, 
2.240, 4.480, 8.960, 17.92, 35.84, 71.68) 



Figure 7. Spectral-index maps of 1803+784 before (left) and after (right) alignment, with the contours of total intensity 
at 1.6 GHz from Fig. 6 superposed. The convolving beam is shown in an upper corner of each image. The shown ranges of 
spectral indices are from —1.50 to 0.50 in both cases. False optically thin emission is visible to the East of the / peak in the top 
spectral-image map. This artefact is absent from the bottom spectral-index map; the "slanting" boundary between the regions 
of optically thick and thin emission near the core seems suspicious, but in fact, the initial direction of the 4.8-GHz jet on small 
scales is to the Northwest (Gabuzda 1999). 
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Figure 8. Cross-correlation functions obtained when the 7.9 & 5.1-GHz / images for 2007+777 (top left) and 2200+420 (top 
right) and 4.8 and 1.6-GHz images for 1803+784 (bottom) were compared. Contour levels are 90, 95, 98, 99.5, and 99.9% of 
the peak values of 0.9927 (top right), 0.9991 (top left) and 0.9185 (bottom). The units plotted along the axes are pixels. 



